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ABSTRACT 

We present high-resolution X-ray spectra from the young supernova remnant Cas A using a 70-ks 
observation taken by the Chandra High Energy Transmission Grating Spectrometer (HETGS). Line 
emission, dominated by Si and S ions, is used for high-resolution spectral analysis of many bright, narrow 
regions of Cas A to examine their kinematics and plasma state. These data allow a 3D reconstruction 
using the unprecedented X-ray kinematic results: we derive unambiguous Doppler shifts for these selected 
regions, with values ranging between —2500 and +4000 kms -1 and the typical velocity error less than 
200 kms -1 . Plasma diagnostics of these regions, derived from line ratios of resolved He-like triplet lines 
and H-likc lines of Si, indicate temperatures largely around 1 keV, which we model as O-rich reverse- 
shocked ejecta. The ionization age also does not vary considerably over these regions of the remnant. The 
gratings analysis was complemented by the non-dispersed spectra from the same dataset, which provided 
information on emission measure and elemental abundances for the selected Cas A regions. The derived 
electron density of X-ray emitting ejecta varies from 20 to 200 cm~ 3 . The measured abundances of Mg, 
Si, S and Ca are consistent with O being the dominant element in the Cas A plasma. With a diameter of 
5', Cas A is the largest source observed with the HETGS to date. We, therefore, describe the technique 
we use and some of the challenges we face in the HETGS data reduction from such an extended, complex 
object. 

Subject headings: radiation mechanisms: thermal — supernova remnants — ISM: individual(Cas A 
(Gil 1.7-2.1) — Techniques: Spectroscopic — X-rays: ISM 



1. INTRODUCTION 

Cassiopeia A (Cas A, Gill. 7-2.1) is the youngest Galac- 
tic supernova remnant (SNR), believed to be a product of 
a SN explosion in ~ 1670 (e.g. Thorstensen et al. 2001; 
Fesen et al. 2006). The SNR is the brightest object in the 
radio band, and is very bright in the X-ray band, providing 
a suitable target for detailed studies across a wide range of 
wavelengths. Cas A shows a wealth of phenomena impor- 
tant for studying the SN progenitor, its explosion mech- 
anism, the early evolution of SNRs, and their impact on 
the interstellar medium. 

The remnant appears as a bright, clumpy ring of emis- 
sion with a diameter of around 3' associated with the SNR 
ejecta, while the fainter emission from the SNR forward 
shock forms a filamentary ring with a diameter of 5'. The 
almost perfectly circular appearance of the SNR is dis- 
rupted by the clear extension in the NE region of the SNR, 
called the jet, and a less obvious one in the SW called the 
counter-jet (e.g Fesen & Gunderson 1996; Hwang et al. 
2000). 

Many of Cas A's bright knots have been identified as 
shocked ejecta, still clearly visible due to the young age of 
the SNR. Observations in the optical band provided the 
first insights into the kinematics and chemical composi- 
tion of these ejecta (e.g. Minkowski 1959; Peimbert & van 
den Bergh 1971; Chevalier & Kirshner 1978). The optical- 
emitting ejecta have been classified into two main groups: 
the fast-moving knots (FMKs), moving with speeds from 
4000 kms -1 up to 15,000 kms -1 , and the slow-moving 
quasi-stationary floculi (QFS), moving with speeds less 
than 300 kms -1 (Kamper & van den Bergh 1976). The 



emission from FMKs, believed to be ejecta, is H-deficient 
(from the lack of Ha emission) and dominated by forbid- 
den O and S emission (e.g. Fesen et al. 2001). This lack 
of hydrogen in FMKs suggests that Cas A was a Type lb 
supernova produced by the core-collapse of a Wolf-Rayet 
star (e.g. Woosley et al. 1993). The emission from QFSs is 
rich in N, which is the reason these knots are believed to 
originate from the circumstellar envelope released by the 
progenitor star and subsequently shocked (or "over-run"). 

While the available kinematic information from the opti- 
cal observations of Cas A is based on observations of thou- 
sands of individual knots, X-ray observations are more dy- 
namically important because X-rays probe a much larger 
fraction of the ejecta mass (more than 4M©) than does the 
optical emission (accounted for less than 0.1M©). Kine- 
matics in the X-ray band have been studied using radial 
velocities (Markcrt et al. 1983; Holt et al. 1994; Hwang et 
al. 2001; Willingale et al. 2002) and proper motion (Vink et 
al. 1998; DeLaney & Rudnick 2003; DeLaney et al. 2004), 
but there is still a need for high-resolution spatial and 
spectral observations that will match the quality of results 
provided in the optical band. Here we present observa- 
tions of Cas A with the Chandra High Energy Transmis- 
sion Grating Spectrometer (HETGS). Cas A is one of the 
rare extended objects viable for grating observations due 
to its bright emission lines, narrow bright filaments and 
small bright clumps that standout well against the contin- 
uum emission. 

2. OBSERVATIONS and data analysis 

Cas A was observed with the HETGS on board the 
Chandra X-ray Observatory in May 2001 as part of 
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the HETG guaranteed time observation program (ObsID 
1046). The exposure was 70 ks, the roll angle was 86° 
West to North and the aim point was at RA = 23 h 23 m 29 s , 
Dec. = +58°48'59"6. The HETGS consists of two grating 
arms with different dispersion directions: 1) the medium- 
energy grating (MEG) arm which covers an energy range 
of 0.4-5.0 keV and has FWHM of 0.023 A, and 2) the 
high-energy grating (HEG) arm which covers 0.9-10.0 keV 
band and has FWHM of 0.012 A (for more details see e.g. 
Canizares et al. 2005). The different dispersion directions 
and wavelength scales of the two arms provide a way to 
resolve potential spectral/spatial confusion problems (sec 
Appendix A; also Dewey 2002). The HETGS is used in 
conjunction with the Advanced CCD Imaging Spectrome- 
ter (ACIS-S). In addition to the dispersed images, a non- 
dispersed zeroth-order image from both gratings occupies 
the S3 chip at the aimpoint, while the dispersed photons 
are distributed across the entire S-array. The zeroth-order 
image, thus, has spatial and spectral resolution provided 
by the ACIS detector. 

The Cas A HETGS event file was re-processed using 
the standard procedures in CIAO 1 software version 3.2.2, 
employing the latest calibration files. Processing of the 
data included removal of the pipeline afterglow correc- 
tion (a significant fraction of the rejected events was from 
the source), correctly assigning ACIS pulse heights to the 
events and filtering the data for energy, status, and grade. 
To retain more valid events, the removal of artifacts (de- 
streaking) from the S4 chip was done by requiring more 
than 5 events in a row in order to destreak it. 

The standard procedures in CIAO for gratings data as- 
sume that the source is point-like. We, therefore, used al- 
ternative software (written in IDL) that basically follows 
the steps of the CIAO threads for gratings spectra, but ac- 
counts for an extended, filament-like source during extrac- 
tion of the PHA spectra and the corresponding spectral 
redistribution matrix files (RMFs). The RMFs are par- 
ticularly important as they relate the photon energy scale 
to the detector dispersion scale of the gratings. We also 
used standard CIAO threads to create auxiliary response 
files (ARFs), which contain the information on telescopes 
effective area and the quantum efficiency as a function of 
energy averaged over time. The resulting ARFs were ex- 
amined for bad columns, and the parts of the spectrum 
where the bad columns are present were ignored in the 
fitting procedure. Details of this "filament analysis" are 
presented in Appendix A. 

For the zeroth-order data, we extracted spectra and 
corresponding ARFs and CCD RMFs with the standard 
CIAO thread ac is spec. For the background spectra we 
tried using the emission-free region on the S3 chip, the 
SNR regions surrounding our bright knots, and the regions 
from the ACIS blank-sky event files. We found no differ- 
ence when using any of these spectra, so we decided to use 
the ACIS blank-sky events. The background spectra are 
extracted with region sizes a factor of 2 larger than the 
source spectra. The zeroth-order model fits were carried 
out binning the data to contain a minimum of 25 counts 
per bin. 

3. RESULTS 



Figure Bl shows the non-dispersed (zeroth-order) image 
of Cas A at the center of each panel and dispersed images 
in different energy bands that contain predominately He- 
and H-like ions, and in the case of Fc also Li- and Be-like 
ions: O+Ne+Fc-L (0.65-1.2 keV), Mg (1.25-1.55 keV), 
Si (1.72- 2.25 keV), S (2.28-2.93 keV), Ar (2.96-3.20 keV), 
Ca (3.75-4.00 keV) and Fe-K (6.30-6.85 keV). The Si-band 
image also indicates the two dispersion axes of the MEG 
and HEG, which are rotated against each other by ~ 10°. 
The Cas A X-ray spectrum is dominated by Si and S lines 
and these bands produce the clearest dispersed images. 
The dispersion angle 9 defines the location of the dis- 
persed photon and varies nearly linearly with wavelength 
A, sinO = mX/p, where m is the order of dispersion and p 
is the period of the grating. Thus, the dispersed images in 
the longer wavelengths (O+Ne+Fe-L band) are spread out 
furthest across the spectroscopy CCD array, while in the 
shorter wavelength (Fe-K band) the images stay tight to 
the zeroth-order. Because of this, the longer-wavelength 
images are more sensitive to velocity gradients and dis- 
tortions are more prominent, as indicated by the larger 
smearing seen in the top two panels. The O+Ne+Fe-L 
image is additionally confused due to the many lines in 
this band (e.g., Fe-L lines and the O Lyman series; see 
Blceker et al. 2001), and a larger relative contribution by 
the continuum. 

Figure B2 shows the regions in the image of Cas A 
selected for this study. These regions have a morphol- 
ogy and brightness that allow the reliable measurement of 
line fluxes and centroids. They are spatially narrow along 
the dispersion direction (north-south) and sufficiently iso- 
lated above the local, extended background emission to 
provide a clear line profile. Figure B3 shows close-ups 
of these regions. For each region we concentrated on the 
four strongest dispersed spectra: the MEG +1 and —1 or- 
ders and the HEG +1 and —1 orders. Because of their 
dominant emission we use the He- and H-like transitions 
of Si and S lines for the analysis of velocity and plasma 
structure in Cas A. For typical SNR plasma densities the 
He-like triplet of the Si and S line shows strong forbidden 
(/) and resonance lines (r) and a comparatively weaker in- 
tercombination line (i). We, therefore, jointly fit our four 
HETG spectra (MEG+1 and HEG+1) with a model con- 
sisting of 4 Gaussians (representing 3 He-like lines and 1 
H-like line for each element) and a constant that accounts 
for the continuum level (for more details see Appendix A). 
Figure B4 shows the four sets of Si spectra from region Rl, 
the brightest of the 17 regions, and from region R9 which 
has the most prominent Si XIV line. Figure B5 shows sim- 
ilarly S spectra from these two regions. In the fit we allow 
only the centroid of the r-linc to be a free parameter, while 
the other three centroids are tied to it, so the relative cen- 
troids are fixed, but the template is allowed to slide. We 
also assume that a single velocity is present in each Cas A 
region, so all 4 Gaussians have frozen narrow widths. In 
other words, if there is no velocity structure in a region, 
the chosen Gaussian width would be the narrowest profile 
for that spatial structure. In addition, the flux ratio of the 
/-line to the z-line (so-called R-ratio) for the Si triplet is 
tied to be 2.63, suitable for the very low density regime of 
the Cas A plasma (e.g. Porquet et al. 2001). Similarly, the 
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/-line flux is tied to be 1.75 times the i-line flux for the S 
triplet (e.g. Pradhan 1982). 

To derive accurate Doppler shifts from measuring line 
centroids we used HEG and MEG spectra binned to 0.01 
A and 0.02 A, respectively. Our Gaussian fits are not sen- 
sitive to measure velocity dispersion as done in the case 
of E0102-072 (Flanagan et al. 2004), where topographi- 
cal changes in the dispersed image between positive and 
negative orders were used to determine intrinsic bulk mo- 
tions. The RMFs we use in the fitting of Cas A HETG 
spectra only include broadening due to spatial structure 
and do not include any further velocity structure degrees 
of freedom in the quantitative fitting. But we can qualita- 
tively distinguish regions with more or less velocity struc- 
ture present. Among the line profiles from our 17 regions 
we find two qualitative types: 1) a narrow profile, where 
the r and / lines are clearly resolved and shifted due to 
the region's bulk motion along the line of sight, and 2) a 
smeared profile, where one of the grating orders appears 
broadened due to high-velocity motions within the region 
relative to its center of mass that causes distortion of the 
dispersed images along the dispersion direction. For each 
region the type of the line profile is indicated in Table Bl; 
Figure B4 shows an example of a double-peak profiles, and 
an example of a clearly smeared spectral profiles is shown 
in Figure B6. 

3.1. Line Dynamics - Doppler Shift Measurement 

Results of Si line measurements and derived Doppler 
shift values are summarized in Table Bl. The spa- 
tial distribution of Doppler velocities is shown in Fig- 
ure B7. Derived values range between —2500 kms -1 to 
+4000 kms -1 . The measured velocity shifts for Si and 
S (not listed here) are similar (see Fig. B8). This is not 
surprising since they arise from the same nucleosynthesis 
layer, and have been found to have the same spatial dis- 
tribution (Hwang et al. 2000; Willingale et al. 2002). The 
uncertainties in derived velocities depend, besides on the 
intrinsic energy resolution of the HETGS, on the number 
of counts detected in the line and the errors associated 
with estimates of the continuum contribution. For the Si 
line the statistical errors for the Doppler shift, based on 
fit confidence limits, range range within 200 kms -1 for all 
regions except for regions R10 and R17 which have errors 
of 650 kms -1 and 360 kms -1 , respectively. 

The HETGS derived velocities of our 17 Cas A regions 
are combined with their spatial location on the sky to 
graphically indicate their 3D location and velocities, as 
shown in Figure B8. Our measured Doppler velocity is 
plotted on the y-axis and the x-axis value is the 2D sky 
radial displacement of the region from the expansion cen- 
ter of Cas A on the sky given by Reed et al. (1995). A 
factor of 0.032" ± 0.002" per kms" 1 is used to relate ve- 
locity to spatial location by minimizing (by eye) the shell 
width needed to enclose most of the data points (dotted 
lines.) The velocity center for the shells of +770 kins" 1 
is taken from Reed et al. (1995). For a distance to Cas A 
of 3.4 kpc (Reed et al. 1995) our factor corresponds to a 
fractional expansion rate of (0.19 ± 0.01) % per year. The 
forward shock location at 153" and the mean reverse shock 
location at 95", as determined by Gotthelf et al. (2001), is 
also indicated. 



3.2. Line Diagnostics - Flux Ratio Measurements 

One advantage of the high-resolution grating spectra of 
our Cas A regions over lower-resolution CCD data is that 
we can investigate the plasma state of individual regions 
using individual emission lines. The He-like K-shell lines, 
like those of Si and S present in Cas A, are the dominant 
ion species for each element over a wide range of tempera- 
ture (e.g. Paerels & Kahn 2003). From the He-like triplet, 
the ratio of the forbidden (/) and resonance (r) lines is 
a useful diagnostic for electron temperature (e.g., the G- 
ratio = (/ + i)/r, Porquet et al. 2001), especially since the 
lines arc from the same ion, which reduces dependence on 
the relative ionization fraction. Lines from different ions 
of the same element are also useful because they eliminate 
the impact of uncertainties in abundance, so e.g. the ratio 
of the H-likc to He-like Si lines in conjuction with the G- 
ratio of the He- like lines can give an accurate measurement 
of the progression of plasma ionization. 

To determine accurate Si line flux ratios from our data, 
we fixed the modelled line locations based on our nominal 
binning and fits described above, and then re-fitted HETG 
spectra using coarser-bin with modified errors. This pro- 
cedure, described in Appendix A, is less sensitive to dif- 
ferences between the shape of the analysis RMF and the 
velocity- modified line shapes. The resulting Si f/r and 
Si XIV/XIII line ratios are listed in Table Bl. We also 
calculated the expected line ratios, employing the non- 
equilibrium ionization collisional plasma model with vari- 
able abundances, VNEI (Borkowski et al. 2001) in XSPEC 
with vneivers version 1.1, which uses updated calculations 
of ionization fractions from Mazzotta et al. (1998). We 
then used these model grids to map our measured line ra- 
tios to equivalent plasma temperature kT e and ionization 
timescalc t = n e t. Figure B9 shows the distribution of 
kT e and r for the 17 Cas A regions. The point for region 
R6 falls outside the displayed ratio range because its Si 
XIV/XIII ratio has an extremely low value. The temper- 
atures range between 0.4 to 5 keV, with the majority of 
regions having a temperature between 0.7 and 1.0 keV; the 
distribution of plasma temperature is shown more clearly 
in Figure B10. The ionization timescale ranges between 
10 10 and 4xl0 n cm -3 s. The measured values for region 
R9 and R12 fall off of the grid in Figure B9 because of 
their extremely high f/r line ratios. To assign kT e and r 
values to region R9, we used the lower error limit which 
falls within the NEI grid. For region R12, we used a value 
within a 20-30% discrepancy of the lower limit. Thus, the 
derived values (e.g., n e , i s hock) for regions R6, R9 and R12 
should be taken with some leeway. 

3.3. Zeroth- order Spectra - Density and Abundance 
Measurements 

In addition to the dispersed data, we also use the non- 
dispersed data (from the central ACIS chip seen in Fig. Bl) 
to obtain information on the abundances and emission 
measure for the individual regions. Obtaining informa- 
tion on global plasma parameters such as electron density 
and elemental abundance requires comparing line to con- 
tinuum in the spectra, and also knowing what are the con- 
tributions to the continuum. The dispersed spectra will 
have contribution from superimposed lines and continua 
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from other image regions, although the superimposed lines 
would not have the right energies (see Appendix A for 
more details). Because of this confusion dispersed data are 
not reliable for measuring line to continuum ratio. There- 
fore, to determine these plasma parameters for each of the 
17 regions, we fitted the non-dispersed zeroth-order spec- 
tra with a single VNEI plasma model (Borkowski et al. 
2001), which allows for varying elemental abundances. In 
these fits we fixed the kT e and r values for each region 
according to the values derived from HETG line ratios. 
Cas A spectra, even of isolated features, are very complex 
and four spectral components have been identified (e.g. 
DeLaney et al. 2004) . These four spectral components are 
rarely present at the same location and most of DeLaney et 
al. (2004) knots show characteristics of a single type, with 
some showing mixed characteristics. To avoid complica- 
tions due to varying column density (-/Vh) and continuum 
levels we ignore the low-energy part of the spectrum and 
consider only the range between 1.1 and 8.0 keV encom- 
passing K-shell lines of Mg, Si, S, Ar, Ca and Fe-K. An 
O-rich plasma is often employed when describing Cas A 
ejecta (e.g. Vink et al. 1999; Hughes et al. 2000; Laming 
& Hwang 2003), since optical observations showed that 
ejecta in Cas A are deficient in H and rich in O and O- 
burning products (e.g. Chevalier & Kirshner 1978). We, 
therefore, assume that O dominates the continuum emis- 
sion, and that O provides many of the electrons instead of 
H and He, as in the typical solar abundance plasma. Thus, 
we fix the O abundance to be a factor of 1000 higher than 
the solar value. The results are, of course, not very sen- 
sitive to the exact O overabundance factor. We also fix 
Nh to 1.5xl0 22 cm -2 , which is found to be an average 
value across the SNR (e.g. Vink et al. 1996; Willingalc et 
al. 2002). In the fit we allowed the normalization and the 
abundances of Mg, Si, S and Ca to vary. The Ar line is 
not included in the VNEI model (vneivcrs version 1.1), so 
we use a Gaussian to model the Ar emission. The promi- 
nent Fe-K XXV line is only present in regions R9 and 
R16, which is not surprising since they are located in the 
area rich with Fe called the " Fe cloud" (Hwang & Laming 
2003). 

We tried different approaches for handling the red-shift 
parameter z in fitting these CCD data. We first froze the z 
values to those derived from the HETGS data, but the re- 
sulting fits were not acceptable, showing a clear mismatch 
between data and model line peaks. These offsets are likely 
ACIS gain uncertainties and other calibration errors in the 
CCD response. Reasonable fits were obtained by manu- 
ally adjusting the z values separately for each data set; the 
values used varied between —0.005 and +0.035. 

Table B2 lists parameters measured from the zeroth- 
order spectral fits of the 17 regions. The first column lists 
the region emission volume Vr, which was taken to be a tri- 
axial ellipsoid whose 2D-projected area corresponds to the 
spectral extraction region shown in Figure B3 with radii 
a and b. For the third axis along the line of sight we take 
the average value of the two observed axes, c = (a + b)/2; 
the volume of the region is then: 

4 

Vr = -7T a b c. (1) 
The second column lists the normalization factor ( "norm" ) 



used in the VNEI model, e.g., 

10~ 14 f 

V norm = J n e n H dV (2) 

Note that the tabulated norm here assumes that the model 
oxygen abundance is set to 1000. The rest of the columns 
list the abundance ratios with respect to oxygen for the 
elements Mg, Si, S and Ca. Using the equations given in 
Appendix B we derive and tabulate in Table B3 some rel- 
evant physical parameters for the regions based on these 
fit results: the region's electron density, the total mass of 
plasma in the region, the time-since-shocked for the region 
ishock = r/n e , in units of years, and finally the fraction of 
oxygen by mass per region. 

4. DISCUSSION 

4.1. Doppler measurements in Cas A 

The measurement of Cas A Doppler shifts in the X-ray 
band was first conducted with the Einstein X-ray Observa- 
tory using the Focal Plane Crystal Spectrometer (FPCS) 
(Markert et al. 1983). The bulk velocities of two regions 
of Cas A, the SE and NW halves, were measured using the 
line centroids of the resolved Si and S triplets. These obser- 
vations found a velocity broadening and asymmetry in the 
X-ray emitting material, with the NW region having more 
red-shifted emission, and the SE region of the remnant 
having more blue-shifted emission. Markert et al. (1983) 
suggested that the asymmetry could be the result of an 
inclined ring-like distribution of Cas A material, possibly 
influenced by the distribution of the mass-loss material of 
the progenitor. Our HETGS spectra of individual Cas A 
filaments reconfirms this global asymmetry trend. The SE 
regions of the SNR appear to be mostly blue-shifted and 
the regions in the NW have the extreme red-shifted val- 
ues. The asymmetry was also found with moderate spec- 
tral resolution of ASCA (Holt et al. 1994). The ASCA 
observations provided a velocity map on the spatial scale 
of V and were derived using the Si line centroids at the 
CCD resolution. 

Doppler velocity maps with a much finer spatial reso- 
lution were produced using Chandra (Hwang et al. 2001) 
and XMM-Newton observations (Willingale et al. 2002). 
Hwang et al. (2001) used the Si XIII resonance (1.865 keV) 
and Si XIV Lya (2.006 keV) line centroids to derive the 
velocity shifts on a 4" spatial scale. While most of the 
SNR showed velocities between —1500 and +1500 kms -1 , 
extreme velocities of —6000 kms -1 were found in the SE 
region, including the region of our Rl, R2, R9 and R13 
for which HETGS data imply velocities between —2500 
and —1000 kms -1 . This discrepancy is not too surprising. 
Although Hwang et al. (2001) argued that they can ig- 
nore the ionization effects on the Si Hea centroid to a rea- 
sonable approximation, they note that high spectral res- 
olution measurements are more desirable to measure the 
energy shifts directly from resolved rather than blended 
lines. Indeed, their Chandra data do not resolve the Si 
XIII triplet, and even the Si XIV and Si XIII Lya lines are 
not fully resolved, introducing some level of uncertainty in 
line centroid measurement. The most recent results from 
the 1 Ms Very Large Project (VLP) Chandra data show 
values more consistent with our HETGS data (Stage et al. 
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2006), but even here absolute gain calibration accuracies 
for the CCD may only be good to 0.5% or 1500 kms -1 . 

XMM Dopplcr maps of Si, S and Fe lines by Willingale 
et al. (2002) show a smaller range in velocity than the 
Si velocity map of Hwang et al. (2001) and, thus, values 
that are at a glance closer to our values in Table Bl. The 
XMM data have been binned to a 20" x 20" spatial grid 
and the spectra were fit with two thermal NEI compo- 
nents representing the ejecta and the shocked component, 
each with a separate energy shift. The Doppler shift val- 
ues derived in this way will depend on the ability of the 
NEI spectral model to predict the line blends combined 
with the uncertainties in the gain calibration of the de- 
tectors. In comparison to our results, it is obvious that 
some of the errors in their Doppler shift values are pro- 
duced because of the spatial averaging over a variety of 
features with very different velocities. For example, in the 
NW region the Dopplcr velocity map of Willingale et al. 
(2002) has a smooth distribution of red-shifted values with 
a 1000 kms -1 or so dispersion, whereas our regions R5 and 
R6 in that part of the SNR show significantly blue-shifted 
values of around —1500 kms -1 , as well as red-shifted ve- 
locities with 1000 kms -1 difference between regions R8 
and R15. Willingale et al. (2002) find that the velocity 
patterns for S are very similar to those for Si; our mea- 
sured Si and S velocities agree with this. The Fe-K ve- 
locities, however, exhibit higher velocities than those of Si 
and S; this has also been found with the 1 Ms VLP Chan- 
dra observations of Cas A (Stage et al. 2006). Note that 
our two regions showing Fe-K emission, R9 and R16, are 
located in the middle (R9) and outside (R16) the range of 
other regions (see Fig. B8). Unfortunately, our HETGS 
data do not have enough counts in the Fe-K band to mea- 
sure shifts in the Fe energy; deeper HETGS observations 
might yield Fe-K velocity measurements with 500 kms -1 
errors, provided the required narrow features are present. 

Probably the best demonstration of the magnitude of 
the difference between the CCD and HETGS measured 
Dopplcr shifts is given by comparing the estimated shock 
expansion rate in Cas A. DcLaney & Rudnick (2003) de- 
rived forward shock expansion measurements of 0.21% 
per year using transverse velocity measurements of Cas A 
knots using Chandra CCD data. DeLaney et al. (2004) 
compared their transverse velocities (3100-3900 kms^ 1 ) 
with that of Willingale et al. (2002) (1000-1500 kms" 1 ) 
and Hwang et al. (2001) (2000-3000 kms" 1 ) obtained 
from Doppler measurements using CCD spectra, and gra- 
ciously ascribe mismatch to possible projection effects of 
the asymmetric remnant. However, our Doppler measure- 
ments imply an ejecta expansion of 0.19% per year, con- 
sistent with the DeLaney & Rudnick (2003) and DcLaney 
et al. (2004) value. This supports their suggestion that 
there might be a dynamic coupling between forward shock 
and ejecta, and they are both part of one homologously ex- 
panding structure. Note that Vink et al. (1998) derived an 
expansion rate of 0.2% per year by comparing the ROSAT 
and Einstein X-ray images from two different epochs. 

The region R17 is the only region outside our reverse- 
shock/forward-shock range and its location at a large ra- 
dius from the expansion center puts it outside the nomi- 
nal forward shock 3D radius as seen in Figure B8. This is 
also a characteristic of the optical FMKs, which are found 



mostly in the NE part of the remnant, with many of them 
located along the jet. Our region R17 is located in the base 
of the jet region (see Fig. B2), and its velocity indicates 
that it indeed might be part of the jet feature. 

4.2. Plasma properties 

Beside the Doppler shifts, plasma temperature and ion- 
ization timescale are the two other plasma parameters de- 
rived directly from our HETGS analysis. Early observa- 
tions of Cas A in X-rays identified two distinct plasma 
components — the cold component with temperature 
around 0.6 keV, and the hot component with tempera- 
ture up to 4 keV (e.g. Becker et al. 1979), and they have 
been confirmed with newer observations (e.g. Vink et al. 
1996; Willingale ct al. 2002). The cooler component is 
associated with the reverse shock traveling through the 
expanding ejecta, and the hotter component is associated 
with the forward shock propagating into the circumstcllar 
material. Plasma temperatures derived from our HETGS 
data have values mostly around 1 keV. Therefore, most 
of our regions have been heated up by a reverse shock 
propagating inward into the supernova ejecta. There is 
no significant pattern to the variation in the temperature 
across the SNR, as shown in Figure B10. The exceptions 
are regions R8, R10 and R17 which do have temperatures 
of 4-5 keV, several times higher than most other regions. 
Regions R8 and R10, which also share similar velocity, 
could, therefore, be associated with circumstcllar material 
and the forward shock. Region R17 could be an X-ray 
counterpart of the optical FMKs found in Cas A jet. Will- 
ingale et al. (2002) derived a map of ionization timescale 
for the cool plasma component in Cas A, which shows some 
variation across the surface of the SNR, but most of the 
SNR has r values larger than 10 11 cm~ 3 s. For our 17 re- 
gions we also find little variation, most of the regions have 
r around a few 10 11 cm~ 3 s. The exception are, again, only 
regions R8 and R10, as shown in Figure B10, which have 
ionization timescale up to an order of magnitude smaller. 

Information on density and abundances for the 17 Cas A 
regions is derived in conjunction with the zeroth-order 
data. The zeroth-order spectra are reasonably fit with a 
single VNEI model with fixed temperature and ionization 
timescale, certainly well enough for our primary purposes 
here to establish relative abundances. An example spec- 
trum is given in Figure Bll for region R9. In the 1.1 and 
8 keV range the spectrum shows a weak Mg lines, strong 
Si, S, Ar and Ca lines, and a weaker Fe-K line, and shows 
reasonable agreement between the data and the model for 
the continuum part of the spectrum. Note that a low- 
energy "up-turn" is also present due to emission from Fe- 
L lines and this part of the spectrum, below 1.1 keV, was 
not used in the fit. The Si XIII Ly/3 line at 2.18 keV is 
generally under-fit in the spectra; this is likely due to sev- 
eral calibration issues each at the 10-20% level (e.g., Ir 
contamination is not included in the HRMA model, there 
are possible zeroth-order HETG calibration errors, CCD 
gain offsets are coupled with steep ARFs.) Future work 
using the higher-count Cas A VLP data set and including 
these calibration effects could usefully extract the He-like 
Ly/3 line flux allowing associated diagnostic ratios (Por- 
quet et al. 2001) and possibly indicating charge-exchange 
processes (Pepino et al. 2004). 
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Electron density for each of the regions, derived as de- 
scribed in Appendix B, is listed in Table B3. The n e ranges 
from around 20 to 200 cm" 3 and it seems to have higher 
values for the blue-shifted regions. A factor of 5 density 
difference has been suggested between the front and the 
back of Cas A (Reed et al. 1995). The time since individual 
Cas A regions have been shocked, given by t s h oc k = i~/n e , 
is also listed in Table B3. These times vary significantly 
from region to region and it seems that the red-shifted re- 
gions have been shocked more recently compared to the 
regions on the front side which have generally larger i s hock 
values (regions R6, R9 and R12 have extreme values that 
should be taken with caution), as shown in Figure B12. 
However, these values could be the result of density dif- 
ferences between the front and the back side of the SNR. 
We do not find any correlation between kT e and i s hock or 
n e and r that would indicate a possible electron-ion equi- 
libration occurring in our 17 regions. 

The other parameters listed in Table B3 are the total 
mass and the fraction of O mass per region. These values 
should be taken with caution, because we assumed in our 
fit that continuum is dominated by O. Note also that our 
total mass estimate does not include Fe or Ar since they 
are not included in our NEI model, so the M to tai values 
could be somewhat modified when these are included. The 
O mass fraction ranges from 0.82 to 0.97, thus confirming 
that, under our assumptions, O is the dominant element 
in these Cas A ejecta features. Cas A is an O-rich SNR, 
as indicated by optical observations, and this leads to the 
further assumption that perhaps there is also a lot of O in 
the regions between the bright Si knots that we see. For 
example, Laming & Hwang (2003) estimate that a den- 
sity enhancement of around 2 coupled with the presence 
of the Si-Ca metals would allow these knots to stand out 
even though surrounded with a pure oxygen plasma. This 
may be another reason that the O-band dispersed image 
appears so smeared out and that we do not detect clear 
O-lines in our HETGS data: there is O emission every- 
where within the shocked ejecta region so we do not see 
discrete O filaments or clumps like those in the Si and S 
band images. 

In order to test the suggestion by Laming & Hwang 
(2003) about the wide-spread presence of O emission, we 
make a few simple estimates. Table B3 shows the mass 
fraction of oxygen to be 0.82 in region Rl. Under the as- 
sumption that the entire SNR volume within a radius of 
100" and 130", for example, is filled with oxygen of the 
density projected for Rl, the total mass of oxygen would 
amount to 85M Q . Earlier studies imply a mass of 4M Q for 
the entire ejecta at most (e.g. Vink et al. 1996). Thus, our 
results seem to either overestimate the oxygen density or 
the density of that particular region is not representative 
for such a large volume. Even if we further assume that 
the density in Rl is enhanced with respect to the ambi- 
ent medium by the factor of 2 to 5 (Reed et al. 1995), we 
would wind up with 40 to 16M Q of oxygen under uniform 
conditions. Introducing a density gradient to the ejecta, 
such as a rise in density towards the center of the remnant, 
gives an increase of at least another factor of 2. Therefore, 
such simple estimates do not seem to support the Laming 
& Hwang (2003) assumption. Alternatively, Willingale et 



al. (2003) derived a filling fraction of 0.009 for the ejecta 
component by assuming a pressure equilibrium between 
hot and cool plasma components, which lead to an elec- 
tron density range of 40-90 cm~ 3 , similar to our values. 
Adding such a factor to our estimate above would then 
allocate much of the metals and oxygen into spaghetti- 
like filaments requiring a total mass of order 1M Q . This 
value represents a more reasonable contribution to a total 
oxygen mass of 2.6M suggested by Vink et al. (1996). 

5. SUMMARY 

High-resolution HETGS observations from the young 
SNR Cas A yield unprecedented kincmatical X-ray re- 
sults for a few bright SNR regions. These observations 
show that high-resolution X-ray spectroscopy is catching 
up with that in optical, IR and UV bands in its ability 
to measure velocities and add a third dimension to the 
data. Unambiguous Doppler shifts are derived for these 
selected regions, with the SE region of the SNR showing 
mostly blue-shifted values reaching up to —2500 kms -1 , 
and the NW side of the SNR having extreme red-shifted 
values with up +4000 kms -1 . This global asymmetry is 
consistent with previous lower spatial or spectral resolu- 
tion X-ray observations. From our Doppler measurements 
we derive ejecta expansion of 0.19% per year, supporting 
suggestion of DeLaney et al. (2004) that there might be 
a dynamic coupling between the forward shock, that is 
expanding at the same rate, and the ejecta. 

Plasma diagnostics using resolved Si He-like triplet lines 
and Si H-like line shows that most of the selected regions 
in Cas A have temperatures around 1 keV, consistent with 
reverse-shocked ejecta. However, two regions, R8 and R10, 
with significantly different temperature of ~ 4 keV, might 
be part of the circumstcllar material. Similarly, the ion- 
ization age does not vary considerably across the remnant, 
except for the two above mentioned regions, which have 
an order of magnitude lower ionization timescale. One 
of the regions, R17, is located outside the forward shock 
boundary and also has a high plasma temperature, both 
suggesting that this region is part of the Cas A NE jet 
feature. 

The analysis of HETGS data was complemented by the 
non-dispersed CCD ACIS spectra from the same obser- 
vation, which allowed us to derive the electron density of 
the X-ray emitting ejecta and the elemental abundances of 
Mg, Si, S and Ar in our 17 regions. The electron density 
varies from 20 to 200 cm~ 3 and does not show any cor- 
relation with ionization timescale. The derived elemental 
abundances of Mg, Si, S and Ca are consistent with O 
being the dominant element in the Cas A plasma. 

We thank John Houck and John Davis for contribu- 
tions in the planning of the Cas A HETGS observa- 
tion, and Glenn Allen, Mike Stage, Kathy Flanagan, 
Tracy DeLaney, Dick Edgar, and Dan Patnaude for use- 
ful discussions. Support for this work was provided by 
NASA through the Smithsonian Astrophysical Observa- 
tory (SAO) contract SV3-73016 to MIT for support of the 
Chandra X-Ray Center and Science Instruments, operated 
by SAO for and on behalf of NASA under contract NAS8- 
03060. 
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APPENDIX 
FILAMENT ANALYSIS 

Filament Analysis Scheme 

The "filament analysis" we used for Cas A HETGS data is very similar to the standard analysis of dispersed grating data 
from a point source, in that it produces a one-dimensional PHA file with associated ARF and RMF for each grating and 
order. The two main additions to the standard processing that adapt it to extended, filament-like sources are described 
briefly below. 

First, adjustments to the event locations are applied to effectively straighten the filament-like source perpendicular to 
dispersion while retaining wavelength accuracy; specifically, the following steps are carried out. The shape of the source in 
zeroth-order is manually traced by a piecewise-linear path that is saved as a set of vertices. The vertex locations are then 
transformed into coordinates aligned with the dispersion and cross-dispersion direction of the particular spectra (HEG or 
MEG) being extracted, see top diagram in Figure B13. Each event, in both the zeroth-order and the dispersed orders, is 
then translated along the grating dispersion direction by an amount equal and opposite to the path offset at that same 
cross-dispersion location. This "shearing" causes the zeroth-order and any dispersed line- images of the feature of interest 
(FOI) to become narrower along the dispersion coordinate, see middle diagram of Figure B13. A PHA file of counts per 
dispersion bin is then created in the usual way by projecting the events along the dispersion axis, see bottom diagram 
Figure B13. Because the dispersed features are effectively narrowed, the ability to detect and resolve discrete lines from 
the feature is improved. 

Second, a companion response matrix, RMF, is created based on the observed, sheared zeroth-order events. For each 
wavelength in the RMF, a subset of events arc selected using the available non-dispersive energy reported by the CCD 
detector and the expected response histogram is created and stored with appropriate offsets in the RMF. Operationally 
this is similar to the effect of the rgsxsrc convolution model available in XSPEC 2 for use with XMM-Newton Reflection 
Grating Spectrometer data. 

This is just one approach to extended source grating analysis (Dewey 2002) and has the advantage of producing familiar 
PHA files which can be analyzed with standard software like ISIS (Houck 2002). Note, however, that it is fundamentally 
an approximation to a fully multi-dimensional, spatial-spectral method and so will necessarily have limitations; some of 
these are implicit in the considerations and techniques described in the following sections. 

Treatment of Spectral Continua 

In addition to line emission from the narrow FOI, there is also continuum emission from the feature, as well as line and 
continuum emission from other parts of the extended source. These give rise to a continuum component in the observed 
PHA spectrum which we discuss and estimate here. 

The observed count rate at a given location on the detector is given by a 3D integral over the source flux as a function of 
wavelength and position on the sky (see eq.(49) in Davis 2001). This integral includes the grating response function, which 
has the property that the location of a dispersed photon includes a continuous dependence on the photon's wavelength. 
Thus, for a given position on the detector there is a set of sky locations and corresponding wavelengths which all contribute 
detected counts in this same location. This is in contrast to the point-source case where the point source introduces a 
delta function in the integral and preserves a one to one mapping of source wavelength to dispersed location in each 
grating-order. The continuum counts seen within a bin in the PHA distribution file will then consist of the sum over the 
spatial regions along the dispersion axis weighted by their flux at the grating-equation-allowed wavelengths. Equivalently, 
the resulting spectrum includes not only the true continuum from the FOI, but the overlapping, shifted spectra from all 
regions along the dispersion axis. 

In practice, the integral may be effectively truncated, e.g., if the spatial extent of the source is moderate. In addition, 
even for very extended sources, the inherent energy resolution of the detector can be used to truncate the integral to a 
limited range in pulse-height, reducing the artificial continuum level. In this case, the observed continuum will be of order 
a factor of R g /RccD greater than the true FOI continuum level, where R g and J?ccd are the effective resolving powers 
of the grating and order-sorting detector, respectively. Note that R g for a given feature varies like 1/E and Rqcd varies 
roughly like \fE, where E is the energy; thus, the R 9 /Rccd ratio varies approximately like E~ 1,5 (or A 15 ). This is a 
small variation, of order ±12% over the 6-7 A range for Si lines. 

As a demonstration of this, a simple MARX simulation was made consisting of a w 4.0" x 10" rectangle emitting at 
1.865 keV embedded in a disk of emission 100" in radius having a uniform spectrum from 1.2 to 2.8 keV. Figure B14 
shows the CCD spectrum extracted from the zeroth-order (e.g., the rectangular region) and the spectrum extracted from 
the MEG dispersed first order. The effective grating resolving power here is R g ~ 75 based on the source full width. 
For the order-sorting effective resolving power, i?ccD, the dispersed extraction was performed with a wide pulse-height 
selection including ±0.15A, giving i?cCD ~ 3.3, based on the energy full width. Together these give R 9 /Rccb ~ 23. In 
the simulation the equivalent width, EW, of the line is 21.4 keV for the zeroth-order case, and 1.29 keV in the dispersed 
case; the continuum level is, therefore, ss 17 times larger in the dispersed data set, which is of the order expected from 
the simple estimate. 

2 http://hcasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/XSmodclRgsxsrc.html 
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Velocity Effects and Fitting 

The features we observe in Cas A show Doppler shifts of up to several thousand kms -1 . For a simple bulk motion of 
the emitting region this just produces an overall Doppler shift to the line wavelength which can be readily measured using 
the filament analysis products. It is also possible that there are velocity variations within the filament itself and this can 
introduce some complications, especially when carrying out standard x 2 -driven fitting. 

In order to study these velocity effects, simple MARX simulations were carried out using a filament source (resembling 
our region Rl) consisting of two sets, upper and lower, of three closely spaced (by precisely 1") parallel line sources. 
The two sets were tilted and offset to produce a wide, kinky filament, shown in panel (a) in Figure B15. In the MEG 
simulations of panels (a) and (b) all the 6 line sources making up the filament all have the same wavelength. These two 
panels show how the shearing of the filament analysis improves the resolution of the dispersed spectra. 

In simulation (c) the upper set of three lines have all been blue-shifted by 1000 kms -1 demonstrating the effect of 
velocity variation along the filament. This causes the dispersed projections to be additionally blurred, having a FWHM 
larger than the (unaffected) zeroth-order. A similar result would arise for the case of turbulent velocity broadening where 
there would be a range of velocities in all parts of the emitting feature. 

The simulation in panel (d) is perhaps the most pathological: here the velocity of the filament varies across the filament, 
along the dispersion direction. The zeroth-order continues to be unaffected by these velocity effects, and now not only 
do the the dispersed order projections differ from the zeroth-order, but they are different from each other as well (see 
HETGS observations of SNR E0102-072, Flanagan et al. 2004). 

The effect of these velocity produced changes in the dispersed line profiles is that the RMF created from the zeroth-order 
is not a good match when there is velocity structure in the feature. For purposes of fitting the line location, the generally 
peaked nature of both the RMF and the dispersed line peaks leads to reasonable centroid fitting and confidence ranges; 
however, the formal x 2 values can be large. 

The mismatch of the zero- velocity line shape used in our RMFs is more problematic when measuring the line fluxes. To 
reduce these mismatch effects we can fix the line locations in the model at the centroid values determined in our nominal 
fitting and then rebin the data and model to a coarser grid. We use binning of order 10 bins, with bin boundaries located 
in the regions between the lines. In this way the total counts in a line region are compared between model and data while 
ignoring line shape differences. Because of the large number of counts in each coarse bin, the errors are now dominated 
by the approximations of our RMFs and systematic calibration errors between the 4 spectra we are jointly fitting, HEG 
m=±l and MEG m=±l. In place of the usual statistical error, we assign a constant error value to the bins of each coarse 
spectrum with a value of 4% of the maximum counts in a bin of that spectrum. This produces fits in which the fluxes of 
the lines are better estimated, as we confirm by examining the fits when data and coarser-fit model are re-plotted to our 
nominal binning, as demonstrated in Figure B16. 



CONVERTING VNEI PARAMETERS TO PHYSICAL UNITS 

This appendix summarizes the conversion of model parameter values into physically meaningful plasma quantities, e.g., 
the electron density, n e , and the mass of each element present, M(Z). A variety of models in the XSPEC library, including 
the VNEI model we use in this work, digest the properties of the emitting plasma into a normalization factor, A norm , and 
a set of relative elemental abundances, Xa(Z). In addition to these, two other source parameters are needed: the source 
distance, 

d [cm] = 3.1 x 10 21 d kpc , (Bl) 

and the volume of the emitting region, 

V R [cm 3 l = I4 s3 (d— — ) 3 . (B2) 

For convenience the conversion from values in kpc, dk pc , and arcseconds 3 , V aS 3, are shown in these equations. 
The actual number fraction of the element Z in the plasma is then given by: 

f(Z) = X A (Z) A modcl (Z) I J2 X A (Z)A modcl (Z), (B3) 

where A mode \(Z) is the reference "solar" number abundance ratio assumed by the model; in our case these are the Anders 
& Grevesse (1989) values. The parameters, Xa{Z), are the usual relative abundance values as input in XSPEC models. 

Because we do not assume that hydrogen dominates the plasma, it is necessary for the electron density to be self- 
consistent with the densities and ionization states of the ions present. The ratio of electron n e to total ion density m is 
then given by: 



& = E^)/( z )> ( B4 ) 



where Q(Z) is the average number of electrons stripped from the ions of element Z, in the range to Z. Ideally Q(Z) 
would be provided by the model; lacking direct access to it, however, we can make a simple approximation to the ionization 
state of the elements in our O-rich plasma by setting: 
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(Z, if Z < 9; 

Q(Z) = \ Z - 2, for 10 < Z < 16; (B5) 
[ Z - 10, if Z > 17. 



The electron density can then be calculated as: 



»e = ^a(^ = 1) (n e n H ) C±) f(z 1 =iy (B6) 

where n# is the hydrogen density and (n e njy) is given by the usual normalization definition 

(n e n H ) = And 2 10 14 X norm /V R . (B7) 

The Xa{Z — 1) factor is included in case the hydrogen model abundance is set to other than 1.0, e.g., to a small value 
like lxl0~ 9 for a pure metal plasma. The density of element Z is then given by 

n 

n(Z) = f(Z) n e I (-£), (B8) 

and other quantities like the mass or mass fraction can be calculated in a straight forward way using the volume, V, and 
appropriate constants (1 amu = 1.66xl0~ 24 g and M©= 2xl0 33 g.) 
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Table Bl 

Measured and derived parameters for 17 Cas A regions from HEG+MEG spectral fits to Si lines: line 

VELOCITIES, PLASMA TEMPERATURE kT , IONIZATION TIMESCALE T, AND Si LINE RATIOS, FORBIDDEN TO RECOMBINATION LINE 

RATIO AND XIV TO XIII LINE RATIO. We ALSO LIST THE TYPE OF LINE PROFILE IN THE LAST COLUMN: " N" INDICATES 
NARROW, NORMAL, DOUBLE PEAKED Si XIII LINE PROFILE PRODUCED BY T AND / LINES, AND "SM" INDICATES SMEARED LINE 
PROFILE DUE TO VELOCITY GRADIENT. NOTE THAT FOR THE THREE REGIONS GIVEN IN PARENTHESIS THE Si RATIOS ARE NOT ON 
THE MODEL GRID (SEE FlG. B9) AND SO THE kT AND T VALUES ASSIGNED ARE APPROXIMATE. 
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4.2 +22 9 xlO n 





.61 + 0.06 


0.12 + 0. 


.03 


n 


R4 


-620 ± 150 


74+ - 11 


2.1±f : |xlO n 





.60 + 0.07 


0.06 + 


.03 


n 


R5 


-1735 ± 118 


77+ - 13 


3.6tf 7 xlO n 





.59 + 0.07 


0.12 + 


.03 


sm 


(R6) 


-1490 ±90 


0.43 (< 1.5) 


3.5xl0 9 '(> 2xl0 9 ) 





.86 + 0.10 


0.001 + 


.027 


sm 


R7 


+3585 ± 135 


1 20+ 20 


4.9+ 2 L3 ,xi0 10 





.41 + 0.04 


0.05 + 


.03 


n 


R8 


+2360 ± 140 




2.2t° : 6xl0 10 





.10+0.04 


0.13 + 


.04 


sm 


(R9) 


-1150+90 


< 1.30 


> 6xlO n 





.70 + 0.17 


0.54+0 


.10 


sm? 


RIO 


+2700 + 650 


4.3ot8i8 


1.3+8- 3 xl0 10 





.11 + 0.01 


0.05 + 


.02 


sm 


Rll 


+310 + 250 


0.99t°-°° 







.54 + 0.08 


0.29 + 


.06 


n 


(R12) 


-850 + 85 


> 0.81 


> 2.6xl0 12 





.97 + 0.13 


0.17+0 


.04 


n 


R13 


-1070+ 140 


0-90±8;?| 


s.et^xio 11 





.54 + 0.08 


0.20 + 


.04 


n 


R14 


+4100+170 


qq+°- 21 


2.4 +2 -ixlO n 





.50 + 0.07 


0.20 + 


.05 


n 


R15 


+760 + 210 




2A +0 1 - 6 4 xl0 11 





.47 + 0.04 


0.17+0 


.03 


sm 


R16 


-1420 + 220 


1 46+ ' 61 


1.2+J : gxlO n 





.37 + 0.08 


0.28 + 


.06 


n 


R17 


-2570 + 360 


4 S0+ 3 ' 70 


7.3+^xlO 10 





.14 + 0.10 


0.96 + 


.18 


sm 



a The velocity values are obtained using full resolution data, in contrast to the other results. 
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Table B2 

Emission volume for 17 Cas A regions, Vr, and measured parameters from VNEI fits to the zeroth-order 

SPECTRA: VNEI NORMALISATION FACTOR, X noIul , AND ELEMENTAL ABUNDANCES WITH RESPECT TO OXYGEN (THE ABUNDANCE 
OF O IS TAKEN TO BE 1000 ITS SOLAR VALUE). NOTE THAT ZERO ABUNDANCE RATIO VALUES DENOTE RATIO VALUES SMALLER 

THAN 0.01. As IN Table Bl, for the three regions given IN parenthesis THE ACCURACY of the derived values 

DEPENDS ON THE MODEL ASSUMED. 



Region 


V R 

(xlO 50 cm 3 ) 


-^norm 

(xl0~ 5 cm" 5 ) 


Mg/O 


Si/O 


s/o 


Ca/O 


Rl 


83 


6 4 +0 - 5 


0.0 


i-s±8:? 


9 fi +0.2 
Z - D -0.1 


1 s +1 - 3 


R2 


44 


n p+0.3 
z - 8 -0.5 


0.0 


o.8l8:i 


i q+0.2 
1 .°-0.2 


Q+ - 8 
u.y_ .4 


R3 


25 


6 fi +0 - 5 

°-°-0.6 


04+ 03 


o.5l8:i 


o.6l8:l 


o.t!8:S 


R4 


28 


7 n+0.5 
' ' U -0.6 


02+ - 02 


o.4i8:l 


o.6l8:l 


o.8l8:2 


R5 


6 


O O+O.l 

d - 8 -0.5 


0.0 


5+ - 1 
u -°-o.i 


0.818;? 


4+ 1 ' 2 

U '*-0.3 


(R6) 


10 


^-0.5 


0.0 


o.3l8:l 


o.6l8:l 


03+ 04 

u - uo -0.03 


R7 


13 


0.918;? 


1+ 01 
u - 1 -o.i 


o.8l8:l 


1 5+ - 1 
1 -°-0.2 


9 9+O.8 
z - z -0.9 


R8 


7 


o.4i8;l 


2+ 01 
u - z -o.i 


o.8l8: 3 2 


ens 


"•^-0.4 


(R9) 


31 




3 +0 1 
u -°-o.i 


7+ - 1 


0.918:? 


5 +0 - 3 

u -°-0.3 


R10 


65 




2+° 1 
u - z -o.i 


1 o +0 - 3 

1 - u -0.2 


i 9+0.4 
1 - z -0.3 


1 +0 - 3 


Rll 


29 


2 9 +0 - 4 

z - a -0.3 


1+ - 4 


o.3l8:l 


o.5l8:l 


1+ - 4 

u - i -0.1 


(R12) 


5 


i-i±°o: 2 6 


0.0 


0-9181 


2 +l fs 


0.0 


R13 


27 


2.8+° ; 5 5 


o.o5l8:8g 


o.5l8:l 


i 1+0.3 
1 - 1 -0.2 


i 1+0.8 
L - L -0.7 


R14 


14 


o c+0.4 
z "°-0.3 


o.o2i8;g 


o.2i8:l 


3 +al 


6 +0 <3 


R15 


16 


9 c+0.1 
z -°-0.4 


0.0 


o.4i8:l 


0-518:? 


0.0 


R16 


109 




o.4i8:i 


o.6l8:l 


o.5l8:l 


o.2l8i 


R17 


87 


0.418;? 




0-618;? 


0-518;? 


0.2I8I 
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Table B3 

Parameters derived from the zeroth order spectra for the 17 defined regions in Figure B2: electron density 

n e , TOTAL EJECTA MASS Mtotal, TIME SINCE REGION WAS SHOCKED t s hock, AND OXYGEN MASS FRACTION M oxygen / Mtotal • As IN 

Table Bl, for the three regions given in parenthesis the accuracy of the derived values depends on the model 

ASSUMED. 



Region 




M to tal 


^-shock 


M 




( cm~ 3 ) 


(xl0- 3 M o ) 


(yr) 


/ Motal 


Rl 


93 


1.3 


71 


0.82 


R2 


81 


0.6 


43 


0.90 


R3 


161 


0.7 


83 


0.94 


R4 


157 


0.7 


42 


0.94 


R5 


235 


0.3 


48 


0.93 


(R6) 


203 


0.3 


0.5 


0.95 


R7 


82 


0.2 


19 


0.88 


R8 


76 


0.1 


9 


0.91 


(R9) 


89 


0.5 


>200 


0.90 


RIO 


32 


0.4 


13 


0.88 


Rll 


98 


0.5 


116 


0.95 


(R12) 


140 


0.1 


>600 


0.92 


R13 


101 


0.5 


112 


0.92 


R14 


131 


0.3 


58 


0.97 


R15 


123 


0.3 


62 


0.96 


R16 


47 


0.9 


80 


0.92 


R17 


23 


0.3 


101 


0.90 




Fig. Bl. — HETG data of Cas A: images of the different line bands. Grating dispersion axes and north direction are marked in the Si-band 
image. Note the smeared out dispersed-order images, especially in the O+Ne+Fe-L band, due to multiple lines and velocity shifts which, 
however, do not affect the central, zeroth-order image. 
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Fig. B2. — Zeroth-order Cas A image from this work. Regions used in the analysis of HETGS Cas A data. These regions are spatially 
narrow and they are isolated sufficiently above the local and extended background to provide a clear line profile for spectral fitting. 
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Fig. B3. — Close up of the regions used in the analysis of HETGS Cas A data. Note that the grayscale dynamic range is not the same for 
all the sub-plots. Some of the regions have been covered with two elliptical regions to follow the HETGS extraction path (sec Appendix A). 
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FlG. B4. — Examples of HETGS spectra of the Si band encompassing four Si lines used in Cas A analysis. The upper four panels show 
the MEGzbl (left) and HEG±1 (right) spectra from region Rl. The model is plotted with the red line. The nominal positions of Si lines are 
marked at the top of each panel. For region Rl the observed lines are shifted to lower wavelength, so this region is an example of a blue-shifted 
region. The lower four panels similarly show the spectra from region R9; this region has a very high Si XIV to Si XIII ratio. Region Rl is 
also an example of a strong, clear line profile produced by a well isolated clump and little of velocity smearing, whereas from region R9 the 
line profiles are messy due to lower contrast between the clump and the surrounding medium and possible velocity smearing. 
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FlG. B5. — Examples of HETGS spectra of S band encompassing four S lines used in Cas A analysis. The upper four panels show the 
MEGzbl (left) and HEGzbl (right) spectra from region Rl. The lower four panels similarly show the spectra from region R9. As in the 
previous plots, the model is plotted with a red line. The nominal positions of S lines are marked at the top of each panel. 
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Wavelength [Angstrom] Wavelength [Angstnom] 



Fig. B6. — Example of velocity-smearing in region R15 due to a velocity gradient in the dispersion direction (see also Fig. B15(d)). Note 
that for both the MEG and HEG spectra the —1 order (upper plots) is narrower and more peaked than the +1 order (lower plots). As in 
previous plots, the model is plotted with a red line. The nominal positions of the Si lines are marked at the top of each panel. 
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Fig. B7. — Doppler velocity values in kms 1 for individual Cas A regions. The south-east SNR regions measured here are all blue-shifted, 
while extreme red-shifts are seen in the north-west of the SNR. 
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Fig. B8. — 3D Location and velocity of the regions. The line-of-sight velocity, as measured by Si lines and S lines, is plotted vs. the 2D 
projected distance on the sky from the nominal expansion center of Reed et al. (1995). The Si velocity values (diamonds) arc given with 90% 
confidence error bars; the velocity values for S (stars) are given just for the best-fit value. An expansion rate of 0.19 % per year relates the 
velocity and distance scales (see text). The inner and outer solid lines indicate the locations of the reverse shock (95") and the forward shock 
(153") in year 2000 and centered on the Reed et al. (1995) velocity center (marked with triangle), V CC ntor = +770 kms -1 . Dotted lines at 
102" and 130" are a guide to show the shell in which our regions generally lie. For reference, the location of the compact central object in 
Cas A (e.g. Chakrabarty et al. 2001) is also shown; its line-of-sight velocity is unknown. 
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Fig. B9. — Si line ratios overlaid onto a NEI grid. Upper panel shows the full range of values for all the regions (region R6 is located off the 
graph), and the lower panel zooms in on the cluster of values in the upper region of the graph. Temperature, kT e , increases along the lines 
which go generally from upper-left to lower right. Ionization time, t, increases along the lines which are nearly horizontal at left and move to 
the upper-right ending at the high-r limiting asymptote. A few kT values are labeled in red, while some r values are labeled in blue. 
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Fig. BIO. — Plot of plasma temperature kT e vs. ionization timescale r = ra e i for 17 Cas A regions. The upper panel shows the full range 
of values, and the lower panel zooms in on the cluster of values. Three regions, R8, R10 and R17, show different properties from the rest. 
Region R6 is represented with two points, representing the mean value (lower left corner in the upper panel) and the upper limits (lower right 
corner in the upper panel). 
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Fig. Bll. — The zeroth-order spectrum from region R9. The spectrum shows a good fit to the Mg, Si S, Ar and Ca lines and to the 
continuum part of the spectrum between 1.2 and 6 keV; Fe lines were not included in the fit and the part of the spectrum below 1.1 keV is 
not used in the fit. 
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Fig. B12. — The time-since-shocked in years is indicated for each region in a similar plot to Fig. B8. The red-shifted regions appear to have 
been shocked more recently compared to the regions on the front side which have longer i a hock values; the low Si XIV region, R6, with 0.5 yr 
is an exception. 
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Fig. B13. — Filament analysis schematic. A path is defined and used to straighten a filament-like feature seen in zeroth-order. A "shearing" 
process is carried out to translate each event location in the dispersion direction by an amount given by its cross-dispersion coordinate and 
the path. An effective RMF and corresponding PHA file are created and can be used in data analysis. 
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Fig. B14. — Elevated continuum level in dispersed spectrum. The zeroth-order and MEG-first order spectra are shown for a simple MARX 
simulation of a monochromatic source embedded in a large region of continuum. The spectra are normalized to unit flux-in-the-line so that 
the equivalent widths, EW, of the lines serve as a measure of the continuum level. The EW for zeroth-order (left, dashed line) is 21.4 keV, 
whereas the dispersed spectrum has a reduced EW of 1.29 keV (right, dotted). Hence, relative to the line flux, the dispersed spectrum has 
an artificially higher continuum level due to overlapping continua from other locations along the dispersion axis. 
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Fig. B15. — Velocity effects on dispersed line shapes. Each of the four panels here is from a simple MARX simulation and filament analysis, 
showing close-ups of the minus order (MEG —1), zcroth-ordcr (Oth), and plus-order (MEG +1) with event scatter-plots and ID histograms. 
The top panel (a) shows the original, simulated filament without shearing applied. In (b) all events have been sheared using a simple path; the 
width of the feature has been narrowed to about 5 pixels wide. In (c) the upper-half of the feature is emitting at a wavelength blue-shifted by 
1000 kins -1 with respect to the lower half. The zeroth-order is unaffected, but the first orders are blurred equally by this velocity variation 
in the cross-dispersion direction. In (d) the source has a velocity variation along the dispersion direction with blue-shift values of: 0, -500, 
and -1000 kms" 1 . In this case the effect is to broaden the minus order (showing the three individual line sources used to simulate the wide 
filament) and narrow the plus order to the point that it is actually narrower than the zeroth-order itself; the zeroth-order is unaffected. 
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FlG. B16. — An example of coarse fitting to determine line fluxes. The top panel shows the data and model fit for the HEG m=— 1 spectrum 
using our nominal binning; the HEG m=+l and MEG m=±l spectra are simultaneously fit as well but not shown. This fitting allows the 
centroids to be accurately measured. In the middle panel the data and model (with line centroids fixed) are binned on a coarser custom 
grid that is aligned with the line locations. Using this coarse binning the spectra are jointly fit with a constant-percentage error term since 
systematic errors now dominate. In the bottom panel the coarse-fit model and data are shown on the nominal grid; note that the modelled 
lines have grown. 



